arXiv:l509.07154vl [hep-ph] 23 Sep 2015 


DESY 15-178 
MITP/15-077 
September 2015 


ISSN 0418-9833 


Inclusive photoproduction of bottom quarks 
for low and medium pt in the 
general-mass variable-flavour-number scheme 

G. Kramer^, and H. Spiesberger^ 

^ II. Institut fiir Theoretische Physik, Universitat Hamburg, 

Luruper Chaussee 149, D-22761 Hamburg, Germany 

^ PRISMA Cluster of Excellence, Institut fiir Physik, 

Johannes Gutenberg-Universitat, 55099 Mainz, Germany, 
and Centre for Theoretical and Mathematical Physics and Department of Physics, 
University of Cape Town, Rondebosch 7700, South Africa 


Abstract 

We present predictions for 6-quark production in photoprodcution and compare with 
experimental data from HERA. Our theoretical predictions are obtained at next-to- 
leading-order in the general-mass variable-flavor-number scheme, an approach which 
takes into account the finite mass of the 6 quarks. We use realistic evolved nonpertur- 
bative fragmentation functions obtained from fits to e'’“e“ data. We find in general 
good agreement of data with both the GM-VFNS and the FENS calculations, while 
the more precise ZEUS data seem to prefer the GM-VFNS predictions. 

PAGS; 12.38.Bx, 12.39.St, 13.85.Ni, 14.40.Nd 


1 



1 Introduction 


The investigation of heavy-quark (bottom or charm) production in ep collisions is a useful 
test of perturbative Quantum Chromodynamics (QCD) since the heavy quark mass provides 
a hard scale that allows to perform calculations within perturbation theory. At leading order 
(LO), boson-gluon fusion, jg —)■ bb is the dominant process for bottom-quark production. 
When the negative squared four-momentum of the initial and hnal electron, Q^, is small, 
the process ep —)■ ebX is treated as photoproduction, in which a quasi-real photon emitted 
by the initial electron interacts with the partons in the proton (direct contribution). For 
~ 0, in addition to the direct contributions, there are also resolved contributions, where 
the incoming photon changes into an initial quark or gluon, which interact with partons 
from the incoming proton. 

Bottom photoproduction has been measured using several different methods by the HI 
and ZEUS collaborations at HERA. In most of the measurements, the cross section 
was obtained using semi-leptonic decays into muons or electrons. All these measurements 
have been compared with predictions based on the fixed-flavor-number scheme (FENS) [^. 
In this scheme, the bottom quark is generated in the hard scattering-process and appears 
only in the hnal state. Predictions for the photoproduction cross section are calculated 
at next-to-leading-order (NLO) taking the hnite mass of the bottom quark explicitly into 
account. Reasonable agreement of the data with FFNS predictions has been obtained. 

In the large-pr region, characterized by pt 3> mh, the so-called massless or zero-mass 
variable havor-number-scheme (ZM-VFNS) [4]-[^ is considered more appropriate. This is 
the conventional parton model approach, where the b quark, considered massless like any 
other incoming parton, is also an incoming parton from the proton and, for the resolved 
contribution, from the incoming photon. This gives rise to additional contributions from 
hard-scattering subprocesses with b quarks in the initial state, besides those with u, d, s 
and c quarks and the gluon (g). Predictions in this approach are reliable only in the large 
Pt- region where terms of the order ml/p^ are negligible. Calculations at NLO automati¬ 
cally resum leading and next-to-leading logarithmic terms oc ln(p^/mj). At the same time, 
all non-logarithmic terms through O(a^) are retained for = 0. With the usual choice of 
renormalization {pr) and factorization scales (/i/ and fip for initial and hnal-state factoriza¬ 
tion, resp.), pr = Pi = pf = rriT = \Jpt + results are dominated by contributions 

from the 6-quark PDF of the proton down to Pt = 0. 

The neglect of the terms of order in the hard scattering cross section is remedied in 

the general-mass variable-flavor-number scheme (GM-VFNS) [^. This theoretical frame¬ 
work combines the FFN and the ZM-VFN schemes. The mass-dependent terms of the 
FFNS are added to the ZM-VFNS by applying subtractions in such a way that in the limit 
Pt ^ oo the hard-scattering cross sections in the ZM-VFNS are recovered. However, with 
the conventional choice of scales indicated above, the results in the GM-VFNS are still 
dominated by the contributions of the 6-quark PDF down to px — 0. Therefore, there is no 
smooth transition from the GM-VFNS at large px values to the FFNS in the small-pT range. 


2 



and the GM-VFNS in its original definition fails to describe the small-pr cross sections. 

The dominance of contribntions with &-qnarks in the initial state at small pr-valnes is 
connected to the fact that this part is treated in the massless scheme. A calcnlation of 
the 6-qnark initiated snbprocesses in a scheme with massive partons is not available for 
photo- and for hadroprodnctioiij^ In addition it wonld be necessary to correctly describe 
the effects of the finite bottom-qnark mass m;, in the DGLAP evolntion. Attempts in this 
direction have been presented in [^. 

The cross section with massless partons is divergent for 0. For a realistic description 

it is therefore necessary to eliminate the contribntion with massless 6-qnarks in the initial 
state in the small-pr region. For inclusive i?-meson production in pp (pp) collisions this 
has been achieved in recent work together with Kniehl and Schienbein by a suitable 
choice of scales pi and fip. In this modified GM-VFNS the FFNS results with the exact 
dependence are recovered at small pt- This way reasonable results were obtained for the 
inclusive S-meson production cross section measured with GDF at the Fermilab Tevatron 
and with LHGb at the GERN LHG. 


It is the aim of this work to apply the same prescription for the choice of the factorization 
scales to 6-quark photoproduction. For the resolved part the modified prescription works 
as described in for hadron-hadron scattering. For the direct part a corresponding 
calculation has not been done yet. We shall do this here and compare our results to recent 
ZEUS and HI cross section data in order to study whether the modifications of the GM- 
VFNS can lead to a better agreement with presently available experimental data. 


Gross sections for inclusive charmed-hadron photoproduction have been calculated in the 
but without the scale modifications described in 


past but without the scale modifications described in M . The corresponding 

experimental data are not available for very small pr, where the modifications would be 
most clearly visible. For somewhat larger px, very good agreement of the various observables 
for inclusive D* photoproduction has been found 13 . 


The outline of this paper is as follows. In Sect. 2 we introduce our strategy for the transition 
to the FFNS and compare our predictions with the ZEUS and HI data. Our conclusions 
are presented in Sect. 3. 


2 Results and comparison with ZEUS and HI data 


In this section, we shall consider a viable framework for inclusive hi, production, where hb 
is any hadron containing a 6-quark. Our aim is to unify prescriptions for the calculation of 
theoretical predictions at small and at large p^-values. We shall compare with the differen¬ 
tial cross section da/dpx measured by the ZEUS and the HI collaborations at HERA. We 
take the 6-quark pole mass to be mb = 4.5 GeV, evaluate at NLO with n/ = 4 

^ For deep inelastic scattering heavy-quark initiated processes at NLO with massive quarks have been 
considered in Ref. [^. 
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and = 328 MeV if iir < rrib and with nj = 5 and = 226 MeV if fiR > mb- For the 


(5) 


MS 


proton PDFs we nse the parametrization CTEQ6.6 
AFG04_BF 


MS 


14 while for the photon PDF we nse 


15 as onr defanlt, bnt we will also consider other PDF choices for comparison. 


We start with resnlts to be compared with the ZEUS data shown in Fig. together 
with onr theoretical predictions. The data have been measnred in varions different analysis. 
Only Eve data points were reported in [^. The other 13 points are from earlier pnblications. 
Their origin can be traced back from Fig. 8 of |^. The ZEUS data are for the differential 
cross section for 6-qnark prodnction as a fnnction of pt (between 5 and 30 GeV) with the 
kinematic constraints < 1 GeV^ for the momentnm transfer sqnared, 0.2 < y < 0.8 
for the inelasticity y = /{xb]S) (where a;Bj is the Bjorken scaling variable and S the 

total energy in the center-of-mass reference frame) and \yb\ < 2.0 for the 6-qnark psendo- 
rapidity. The data points at px = 9.25, 12.5, 17.5, 22.5, 28.75 GeV from have in average 
smaller experimental errors than the earlier measnrements since they are measnred with the 
largest integrated Inminosity of 133 pb“^. Note that not all points are consistent with each 
other, i.e. for some points the systematic errors are presnmably larger than given in [^. 
In most of the previons measnrements of bottom photoprodnction at HERA, the cross 
section was determined nsing semi-leptonic decays into mnons or electrons. In the b- 
qnark photoprodnction cross sections were determined for the production of jets containing 
5-quarks. The cross sections for the production of 6-jets have been converted into cross 
sections for 6-quark production with the help of the FMNR program [^. In Fig. 8 of 
Ref. the data have been compared with predictions of the FFNS (^. 

In Fig. (upper left frame), we compare the ZEUS data with NLO predictions in the GM- 
VFNS. Our framework has been described in detail in Refs. 01 for photoprodnction 
of D* mesons. The full line in Fig. (upper left frame) shows the result for the old 
default choice of scales, pi = ^irriT with = 1 for i = R, J, F where R, /, and F denote 
renormalization, initial-state and final-state factorization scales. The dashed lines represent 
an estimate of the theoretical error by varying up and down by a factor of 2. 13 data 
points agree with the theoretical prediction within experimental and theoretical errors, 5 
points lie outside this range. 

The central curve in Fig. (upper left frame), as well as the curve for the upper error 
estimate, show the characteristic behavior of a monotonic increase of da/dpx in the limit 
Pt 0. This behavior is caused by the scale choice with > 1. The curve for the 
lower error estimate with = 0.5, however, shows a turn-over towards small px with a 
maximum at pt — 3 GeV. This is the same characteristic behavior as in the FFNS 
and is caused in our case by the fact that at around px = 7-7 GeV the factorization scales 
pi = = mx/^ fall below mb where the 6-quark PDF and FF are zero. The same 

typical turn-over was also found in our previous study of inclusive R-meson production for 


pp and pp collisions 10 


Actually, the 6-quark cross section is calculated from the R-meson cross section using 
an evolved fragmentation function (FF) for 6 —R taken from 19 and dividing by the 
branching fraction Rr(6 —)■ R) = 0.40. The effect of the evolved FF is visible at larger 
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Figure 1: Differential cross section da/dpT for bottom-quark production as a function of 
the transverse momentum pt of the bottom quark compared with ZEUS data [^. The 
upper left frame is for the original GM-VFNS choice of default scales Pr = pi = Pf = 
the upper right frame shows results with the new default scale pr = mr, pi = Pf = mT/2. 
For comparison we also show in the lower panel results of the FFNS with pji = pi = rriT 
(no FF). 
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values of pr > 20 GeV where the cross section da/dpr is smaller than without a FF. We 
emphasize that a consistent theoretical prescription of the GM-VFNS requires to include 
a FF describing the 6 —)■ hj, transition since final state singularities have to be subtracted 
from the bare cross section in analogy to the calculation in the ZM-VFNS. 


Unfortunately all the ZEUS measurements are for Pt > 5 GeV, so that any turn-over of the 
cross section da/dpT towards = 0 is not measured. All data points agree reasonably well 
also with the GM-VFNS predictions with the original default scale Pr = = Pf = 


A monotonic increase of the cross section da/dpx for the production of a massive quark 
towards small Pt is unphysicah From our previous work on i?-meson production in pp and 
pp collisions we know how the choice of the factorization scales must be modihed in 
order to obtain cross sections for pt’ —)■ 0 which show the expected behavior. Therefore, 
we choose as the new default scales pr = and pj = pp = rriT/^ such that 


as m 
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both the 6-quark PDF and the FF for the b ^ transition vanish for pp < already 
at non-zero pp- With ^r = = ^p = 0.5 as the new default scale choice, we estimate 

the theoretical error in the usual way by varying the renormalization scale parameter pR 
by a factor 2 up and down about the default scale. We have not introduced an extra 
prescription to freeze pr below mb because the choice of pr is not related to switching off 
6-quark initiated subprocesses. The resulting error band is shown in Fig. (upper right 
frame) together with the full curve for the new default choice and the ZEUS data [^. The 
agreement with the data is only marginally better than in the upper left frame of Fig. [T] 
Since the data have pp values larger than 5 GeV the turn-over of da/dpp towards pr —t 0 
is not tested by the ZEUS measurements. 


In the lower part of Fig. [T]we compare the ZEUS data with results obtained in the FFNS. 
Here we have hxed the scales to Pr = pi = mp. The FFNS prediction is only slightly below 
the GM-VFNS calculation at low pp, but lies above it by about a factor of two at the largest 
pp. The difference at large pp is mainly due to the fact that the FFNS calculation does not 
include an evolved FF. We observe that the FFNS calculation over-shoot the high-precision 
data points of ZEUS. 


Similar data are available from the HI collaboration at HERA. The most recent publication 
describes results for photoproduction of 6-quarks measured in the process ep —)■ ebbX. 
The decay channel bh —)■ eeX is selected by identifying the semi-electronic decays of the 6- 
quarks. The production cross section is measured in the kinematic range where the photon 
virtuality is small, < 1 GeV^, the inelasticity is restricted to the range 0.05 < y < 0.65 
and the pseudorapidity of the 6-quarks is in the range \yb\ < 2. The measured cross sections 
are converted into single-inclusive 6-quark production cross sections da/dpp for four pp bins 
in the range 0 < pr < 30 GeV. We have calculated the ep —)■ ebX cross section da/dpp 
as a function of pp in this pp range and with the same kinematic constraints as in in 
the original GM-VFNS with default scales ^r = ii = ip = ^ and with the new default 
scales ^R = = ^F = 0.5. The results, together with scale variations as above to obtain 

an error estimate, are shown in Fig. together with the HI data. The upper left panel 
of Fig. 1^ shows results for the original scale setting, the upper right frame is for the new 
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Figure 2: Differential cross section da/dpT for bottom photoproduction as a function of the 
transverse momentum pt compared to HI data [^. The upper left frame is for the original 
GM-VFNS choice of scale parameters Pr = = Pf = ttit and the upper right frame for 

the new default scales pr = tut, pi = Pf = rriT/^. FFNS results with Pr = pi = rriT (no 
FF) are shown in the lower frame. 
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scale choice. The experimental data for pt > 5 GeV agree with the predictions in both 
schemes. Only the cross section da/dpT in the lowest-p^ bin, 0 < < 4.65 GeV, is lower 

than the predictions in both cases. In this p'r-bin the measured cross section is smaller than 
the theory prediction by approximately a factor of 6 in the upper left frame of Fig. [^and 
by a factor 3 to 4 in the upper right frame of Fig. Although the pr-dependence of the 
differential cross section da/dpr exhibits the expected turn-over at pt — 2 GeV if the new 
default scales are chosen, the predicted cross section is still larger than the experimental 
data, even when taking into account the experimental uncertainty. 

A comparison of the data with the FFNS prediction is shown in the lower frame of Fig. 

again choosing scales as Pr = Pi = rriT and without folding with a FF. As above, the 
FFNS results are slightly below those of the GM-VFNS at low px, but above at large px- 
We note that the NLO corrections in the FFNS are positive and lead to an increase of 
the cross section by about a factor of 1.5 at low px- The HI data with its relatively large 
errors are unable to discriminate between the two schemes. We should remark that the 
HI data are presented as a function of the average transverse momentum of two produced 
bottom quarks, while our calculation is for the one-particle inclusive cross section where 
the independent kinematic variable is the transverse momentum of one observed 6-quark. 
This may be the reason that the theoretical predictions shown in do not agree with our 
FFNS evaluation. 


In order to study the influence of the choice of the proton PDF we present also results 
for two other PDFs: GTIO 20 and HERAPDF1.5 (NLO) 1^. The results are shown in 


Fig. (left frame) as ratios to the prediction obtained for the default choice GTEQ6.6. 
The differences to GTEQ6.6 are very small, in particular for GTIO, as to be expected. For 
HERAPDF1.5 the difference is somewhat larger. In the range 2 < p^ < 20 GeV the cross 
section with HERAPDF1.5 is larger than for GTEQ6.6 and for larger px values, px > 20 
GeV, the ratio is smaller than 1 and decreases towards 0.9 for px = 30 GeV. 


We include also a comparison with two different choices of the photon PDFs, shown in the 
right frame of Fig. [3 Again, we present the results as a ratio of cross sections and choose 
the GRV HO p)6,17 photon PDF as an alternative to the AFG04 parametrization. The 


cross section ratio is very close to 1. Only for the smallest considered px value the ratio 
reaches values slightly above 1.06. 


As is well-known, a complete calculation of the photoproduction cross section has to take 
into account two parts, the direct contribution and the resolved contribution. It is inter¬ 
esting to study whether both parts show the same dependence on px- To see this in some 
detail we have calculated these two parts separately, denoted / dpx and da'^^^/dpx- 
Results for the HI kinematic conditions are shown in Fig. For the sake of a more 
clear visibility we display the ratios of da^/dpx and da^°^ / dpx where x = dir, res and 
da^°^/dpx = /dpx + da'^^^/dpx is the sum of the direct and resolved contributions as 
shown Fig. (upper right frame). From Fig. |^we see that the resolved cross section varies 
between 6% and 15% for px < 5 GeV, while it stays near 6% for the larger px values. 
This means that even when the resolved part would vanish for all px the direct part would 
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Figure 3: Left frame: ratios of the differential cross sections da/dpT for bottom pho¬ 
toproduction using different proton PDF parametrizations with respect to the default 

the dashed curve for HERAPDF1.5 1^ 


CTEQ6.6 14 PDF. The full curve is for CTIO 
Right frame: ratio of the differential cross section da/dpT for bottom photoproduction as 
a function the transverse momentum px for the photon PDF GRV HO 16 ,^ with respect 
to the photon PDF AFG04 15 . In both hgures, the kinematic constraints of the ZEUS 


data and the new default scales pr = h/ = = vnxj^l have been used. 


change by less than 15%. This is too small to be identihed unambiguously in the data. In 
particular, a different px dependence of the direct part can not be made responsible for the 
low cross section measured by HI at = 2 GeV. 


3 Conclusions 

We have presented predictions for 6-quark photoproduction at HERA at next-to-leading 
order in the general-mass variable-flavor-number scheme. In contrast to previous calcula¬ 
tions, we have hxed the factorization scale parameters in such a way that 6-quark initiated 
contributions are eliminated at low transverse momenta. With this new prescription we 
revocer the typical low-px behaviour of a hxed-flavor-number scheme. The comparison with 
experimental data from the ZEUS collaboration shows reasonable agreement over the whole 
range of transverse momenta, while we observe a discrepancy in the lowest-pT^ bin of the 
HI measurement. 
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Figure 4: Ratios of da^^'^/dpT and da'^^^/dpT to the total cross section da^°^/dpT as a 
function of the transverse momentum px of the 6-quark. 
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